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Abstract 

A procedure is proposed which accounts for final state interaction correc- 
tions for near threshold meson production in nucleon-nucleon scattering. In 
analogy with the Watson-Migdal approximation, it is shown that in the limit 
of extremely strong final state effects, the amplitude factorizes into a primary 
production amplitude and an elastic scattering amplitude describing a 3 — > 3 
transition. This amplitude determines the energy dependence of the reaction 
cross section near the reaction threshold almost solely. The approximation pro- 
posed satisfies the Fermi- Watson theorem and the coherence formalism. Ap- 
plication of this procedure to meson production in nucleon-nucleon scattering 
shows that, while corrections due to the meson- nucleon interaction are small for 
s-wave pion production, they are crucial for reproducing the energy dependence 
of the rj production cross section. 
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1 Introduction 



The cross sections for 7r° and r) meson production via the pp — > ppir and pp — > ppr? 
reactions, at energies near their respective thresholds, exhibit a pronounced energy 
dependence which deviates strongly from phase space |], §]. Such deviations, which 
occur in other hadronic collisions as well0, are most certainly due to final state 
interactions(.FS7) between the reaction products. While FSI corrections due to 
the short range nucleon-nucleon(iViV) force between the outgoing protons seem ade- 
quate to reproduce the observed energy dependence |I[ and the cross section scale [Q 
for pion production, they fail to do so for the rj production^. Particularly, several 
model calculations |5], ||, 0] have considered the NN — > NNrj through a perturbative 
approach but none of them reproduced the energy dependence near threshold, al- 
though FSI corrections were made to account for the proton-proton interaction. It 
is the purpose of the present work to show that while the meson-nucleon interaction 
influences s-wave pion production only slightly, its effects are crucial for reproducing 
the energy dependence of the pp — > ppt] reaction. 

There exists no theory of final state effects in the presence of three strongly inter- 
acting particles and our first objective will be to develop a procedure that accounts 
for both the NN and meson- N forces. To accomplish this task we consider a pro- 
duction process, ab — > cde. A full dynamical theory of such a process would require 
the solution of two and three-body scattering problems. Particularly, in the Faddeev 
formalism || the full transition amplitude is decomposed as a sum of three terms de- 
pending on which of the three final particle pairs interacts last. The evaluation of 
these terms, the so called Faddeev amplitudes, requires the solution of a set of coupled 
integral equations. Such a procedure should resolve the problem at hand but, in view 
of the scarcity of knowledge about the rjN interaction it may turn to be unreliable 
and rather long and tedious to apply. In what follows we develop an approximation 
along exactly the lines of the Watson's approach || for two body reactions, i.e., we 
look for a separation of the energy dependence due to FSI from those of the primary 
production amplitude. Such an approximation may prove useful and easy to apply in 
analyzing near threshold meson production data from NN scattering. 

According to the Watson-Migdal FSI theory for two body processes 0, the energy 
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dependence of the full transition amplitude is governed by rescattering of final state 
particles. The full transition amplitude is approximated by, 



where the primary production amplitude, M\f , is assumed to be a smooth and slowly 
varying function of energy and t^j is taken as the free two-body elastic scattering 
amplitude in the exit channel. Thus the energy dependence of the full matrix element 
is determined mostly by t ff. A decomposition as such of the amplitude into a primary 
production amplitude and a subsequent FSI correction term is meaningful only if 
iff is sufficiently large so that FSI dominate the reaction amplitude. This requires 
that the distortion of the final state wave function (i.e. the deviation from a simple 
product of free particle wave functions) or alternatively, the sticking factor (i.e. the 
probability of the particles to find each other's vicinity) in itself be proportional to 
t ft . Albeit, the particles spend a great deal of time close together, as should be the 
case in particle production processes close to threshold. Otherwise, the distortion 
factor of the final state (and obviously the full matrix element) must be calculated 
by using more precise methods 0. 

The complete transition amplitude of a production process can be decomposed in 
the form shown in fig. 1, where the first term (the block M^ n ') represents the primary 
production amplitude which includes all possible inelastic transitions contributing to 
the process, while all other terms represent the contributions from rescatterings in the 
entrance and exit channels (the I SI and FSI blocks). The third term, for example, 
corresponds to a process in which the primary production takes place as if there were 
no FSI and only subsequently is distorted by the Coulomb and nuclear short range 
interactions between the reaction products. At energies near the reaction threshold, 
where the relative energies of particle pairs in the final state are sufficiently small 
and FSI corrections are dominating, such a term is expected to be more important 
than others. In analogy to the Watson-Migdal approximation we make the ansatz 
that in the limit of weak ISI and strong FSI , the transition amplitude of a three- 
body reaction can be approximated by an expression similar to eqn. 1, with the FSI 
correction term being replaced by the on mass-shell elastic scattering amplitude of a 
3^3 (three particles in to three particles out) transition. This amplitude is denoted 
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by t!$P and, can be evaluated using the Faddeev formalism! 10] or others such as the 



Weinberg method|jrTJ. Certainly, it is not a measurable quantity as such a 3 — > 3 
transition is not easy to realize experimentally. It is shown in sect. 2 though, that 
at energies near particle production threshold, can be estimated from two-body 
scattering data without having to solve the full three-body dynamics. 

Within this approximation, the transition amplitude of a three-body reaction is 
coherent in terms of the interactions between the various final particle pairs, just 



as one anticipates based on general quantum mechanical considerations |T2|. Thus 
the effects due to a given pair's interaction is distributed over the entire transition 
amplitude and may influence strongly the energy dependence of the cross section. In 
sect. 3 we apply the procedure developed to production of mesons in NN scattering 
and show that although the rj — N interaction is rather weak with respect to the NN 
interaction it makes remarkably strong contributions. Discussion of the results and 
conclusions are given in sect. 4. 

2 Theoretical Perspective 
2.1 Two-body reactions 

Prior to extending the Watson-Migdal approximation to three-body reactions let us 
recall first the more familiar problem of FSI corrections in two-body reactions. As- 
sume that the full interaction between particles separates into v — w + s, where w 
and s stand for weak and strong terms respectively. Here w is the primary interaction, 
such that if it were zero, the process in question would not occur. The remaining part 
of the interaction s is responsible for rescatterings in the entrance and exit channels. 
We further assume that s is the only strong channel energetically allowed. 

We would like to find how FSI affects the transition amplitude in a process in 
which the primary interaction is relatively weak. The full transition amplitude in an 
appropriate eigenchannel is defined as, 

Tit = <*$l |*«> , (2) 

where ^efi(/) are two-particle scattering wave functions. They correspond to solutions 
of the Schrodinger equation with a Hamiltonian that contains strong s-interactions 



only and satisfy the boundary conditions that ty^j (^iij) tends asymptotically to 
free two-particle wave function (p oi (0 O /), plus outgoing (ingoing) spherical waves. 



These can be written in the form[|10 



= lim [ 1 + G 02 (E ± ie) fg>] OA ,X — i, f . (3) 

Here Gq2(E ± ie) is a free two-body Green's functions, and tyf and vf l f are two- 
body elastic scattering amplitude for the entrance and exit channels, respectively. 
Substituting eqn. 3 in 2 leads to the full transition amplitude, 

f if = tiff* + t^G^M^f + AffiGafiff + tH l) G 02 M^G 02 i^ . (4) 

Note that this expression is exact and contains both ISI and FSI . In the Watson- 
Migdal approximation, 

Tit « I M( 0o/ I Mif |^k>< o ^g) ' (5) 

where the two overlapping integrals {4>oi\^a}) an< i (^eij I 00/) are on mass-shell S 
matrix elements which contain ISI and FSI corrections, respectively. To cast the 
amplitude, eqn. 5, in the form given by the Watson's approximation, eqn. 1, one 
has to replace by <f) 0i so that Sa is equal to unity and, in the limit of extremely 
strong FSI , take the on mass-shell instead of the Sff matrix elements. 

It would be instructive to repeat these arguments using a diagrammatic language. 
We associate the weak primary interaction w with the diagram 2a, and the strong s 
interaction with the diagram 2b of figure 2. As usual the disconnected diagram 2c 
describes noninteracting particles. It is evident that the amplitude Tif is the sum of all 
diagrams that can be constructed from the elementary diagrams 2a and 26. To single 
out the ISI and FSI blocks and obtain the full matrix element in a diagrammatic 
representation similar to the one of fig. 1, we first sum in the initial and final channels 
all diagrams that contain strong interactions only (see fig. 3). These sums yield 
the elastic amplitudes (A = z, /) which we identify with the appropriate ISI 
and FSI blocks, respectively. The sum of all diagrams that start or end with a 
weak interaction (diagram 2a) forms the block which describes the primary 

production process. When rescatterings in both of the initial and final states give 
sufficiently large amplitudes, the Migdal- Watson ideology tells us that the 4th term 



on the rhs of eqn. 4 is dominating. When FSI (I SI) effects are small, then the second 
(third) term dominates. Thus the ISI and FSI effects can be separated from the 
full transition amplitude by summing diagrams in a specific order and subsequently 
derive the Watson-Migdal approximation by selecting an appropriate dominant term 
in the limit of extremely weak ISI and extremely strong FSI . We apply this same 
reasoning to obtain the Watson-Migdal approximation for three-body reactions. 

2.2 Three-body reactions 

Let T 23 be the transition amplitude of a production process in which an initial two- 
body process goes into a final three-body state. For conciseness assume all particles 
to be distinguishable and spinless. This restriction simplifies the mathematics but 
can be easily removed to treat more general cases. Furthermore, assume that each 
pair in the three-body final state has a given angular momentum so that a unique 
elastic scattering state is assigned to each pair simultaneously. This is very often 
the case at energies near threshold where, only one partial wave dominates for each 
pair namely, A = 0. Following the discussion given above for two-body reactions 
the interactions are divided into two categories: (i) elastic interactions leading to 
rearrangement of particle (inner) quantum numbers but not to the production of any 
additional hadrons, (ii) inelastic interactions which are responsible for the production 
(or annihilation) of particles so that particle numbers in the initial and final states are 
different. Likewise, we call elastic those diagrams which involve elastic interactions 
only and having equal number of legs to the left hand and to the right hand sides. 
Those which begin and end with inelastic interactions and having different number 
of legs to the left and right hand sides are called inelastic. 

The transition amplitude T 2 3 of the reaction under considerations is the sum of 
all diagrams which start with two legs on the left and end with three legs on the 
right, that can be assembled using all possible elastic and inelastic diagrams. To 
decompose T23 into four terms as depicted in fig. 1, these diagrams need be organized 
in a specific order before suming them. First note that by summing all two-particle 
elastic diagrams in the entrance channel and all three-particle elastic diagrams in 
the exit channel, one obtains the 2 — ► 2 and 3^3 elastic transition amplitude of 
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fig. 3. These elastic amplitudes which are denoted by t%P and T33 , contain only 
elastic interactions and as such, they are only parts of the complete 2^2 and 
3 — > 3 amplitudes. Although at energies above particle production threshold, two 
and three-particle channels may well be opened in both the initial and final states, 
these two amplitudes are calculated by taking into account only elastic interactions. 
All inelastic and quasi-elastic processes are confined into M^ n \ Note also, that the 
2^2 elastic block contains connected diagrams only but, the 3^3 block contains 
disconnected diagrams as well (see fig. 3). 

The set of all diagrams contributing to the entire transition amplitude T23, is now 
separated into four partial sums each corresponding to one of the terms in the block 
diagram of fig. 1 : (i) the block is the sum of all inelastic diagrams only, (ii) the 

ISI -M^ term is the sum of all diagrams which start with a 2 — > 2 elastic diagram 
on the left and end with an inelastic diagram on the right, (iii) the M^-FSI term 
is the sum of all diagrams which start with an inelastic diagram on the left and end 
with a 3 — > 3 elastic diagram to the right and, (iv) the ISI -M^ in '-FSI term is the 
sum of all diagrams which start with a 2 — ► 2 elastic diagram on the left and end with 
a 3 —>■ 3 elastic diagram on the right. Obviously, these four partial sums exhaust the 
entire set of diagrams that may contribute to the complete transition amplitude. 

This same presentation of T23 as the sum of the terms of fig. 1 can as well be 
derived formally. By definition, the transition amplitude is, 

T23 = (^ll\M {m) \^ e tl) , (6) 

where as defined in eqn. 3, stands for a two-body elastic scattering wave function 
and ty^l is a three-particle elastic scattering wave function corresponding to the exit 
channel. It is a solution of a three-body Schrodinger equation with a Hamiltonian 
having only two-body and three-body elastic interactions, and satisfies the boundary 
conditions that ^3 tends asymptotically to a free three-particle wave function, 003, 
plus outgoing spherical waves. In terms of the elastic scattering amplitudes defined 
above one has the expression ]T0|], 



= lim [ 1 + G 03 (E + ie) fjf ] O3 • (7) 
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Using this with eqn. 3 in eqn. 6 implies that, 

f 23 = M (m) + fipGnM™ + M (m) G 3T 3 ( 3 e0 + t^G^M^G^f^ . (8) 

Formally, this expression provides an exact solution of the problem and has a structure 
identical in form to that given in eqn. 4, for two-body reactions. A notable feature of 
this expression is that all inelastic transitions are confined in the block M^ n \ while 
all possible elastic rescatterings in the entrance and exit channels are contained in 
t\i2 and T33 \ 

Let us now consider FSI corrections a la Migdal — Watson. In the limit of 
extremely weak I SI one may replace by 002 in eqn. 6. This makes the second 
and fourth terms on the rhs of eqn. 8 disappear. Furthermore, in the limit of very 
strong FSI the third term on the rhs of eqn. 8 dominates and we may approximate 
the transition amplitude by, 

T 23 ^(0o2|M (m) |0o3)T 3 { 3 e/) . (9) 

We may argue that the matrix element (4>o2\M ( - m ^\(f) 03 ) is a smooth and slowly varying 
function of the energy and momenta and hence the energy dependence of the ampli- 
tude, eqn. 9, is determined almost solely by the elastic amplitude T^ l \ Although 
in form eqn. 9 resembles the Watson-Migdal approximation for two-body reactions, 
eqn. 1, there exists an essential difference between the two cases which makes eqn. 
9 more difficult to apply. This is because the amplitude T^, contrary to tffi , is 
practically not a measurable quantity. 

The two-body initial state, ^^ 2) can be calculated using distorted wave approach 
and the effects of slowly varying nonvanishing I SI can be included by replacing 
in eqn. 9 with the matrix element {(p^M^l^!^^) . With this modification in mind 
eqn. 9 is more suitable to apply to data. Certainly, ISI must be included in a full 
dynamical description of the process as they influence the cross section scale through 
the matrix element M^ n \ 

2.3 The amplitude for 3 -> 3 

In the Faddeev formalism is written as a sum of three Faddeev amplitudes, 

T 3 { 3° = ET J • (10) 



These satisfy the set of coupled integral equations |TIJ, 

oo 

T = t + FT =J2 ( FF ) n [t + Ft 
where we have used the notation, 



(11) 
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(12) 



and final particles are labeled 1,2 and 3. Here t,- are two-body scattering amplitudes 
in the three-body space. In momentum space and adopt the convention that (j, I, k) 
are always cyclic, 



(PjiPhPk I tj | p',pj,p' fe ) = (pi,p k I - tj) I Pi,Pfe) 5(Pj - p'j) 



(13) 



with tj being a solution of the ordinary two-body Lippmann-Schwinger equation of the 
energy variable E — €j. E is the total three-body energy and 6j the energy appropriate 
to momentum pj. From the definitions in eqn. 12, 



T'=tj + (T l + T k ) Goati , J^l^k^j , i = 1 - 3 



(14) 



As pointed by Amado ||12| , the contributions from the various pair interactions add 
coherently and the effects of a given pair's interactions is distributed over the entire 
amplitude. From eqn. 14, T J is the sum of all terms contributing to T 3 g which end 
with tj but, it does not include all contributions of the I — k interaction. Note that 
the first term on the rhs represents the I — k two-body scattering amplitude with the 
jth particle being a spectator and as such it corresponds to a disconnected diagram. 
In the center of mass of the / — k pair and for a given partial wave A, tj has the elastic 
scattering phase 5\(j). The second term in eqn. 14 represents the contributions 
from two or more rescatterings and involves completely connected diagrams only. 
These play an important role in the coherence formalism [|i~2]j through interference 
with similar terms involving the j — I and j — k pair interactions. Therefore, in 
order to account for the effects of any pair's interaction the entire amplitude must be 
constructed. Adding the sum (T l + T k ) to T- 7 of eqn. 14, one obtains three equivalent 
forms of the entire amplitude, 



T 3 { 3 e0 =tj + (T l + T k ) [1 + G m (E) tj] , j^l^k^j , i = 1 - 3 



(15) 



In each of the forms in eqn. 15, there appears the factor (1 + G 3 tj) with the tj being 
half on-shell. At a given I — k relative momentum (q^) in a given I — k partial wave 
(A) and following the arguments in ref.[Q, these half-on-shell amplitudes, the factor 
(1 + G 3 tj) and, the entire amplitude have the phase 6\(]). For this to happen, say in 
a region where the expression (T l + T k ) G 3 tj may vary rapidly with energy it must 
have a part which cancels the variation of (T l + T k ) so as to give to their sum the 
behaviour of tj. This general observation on three-body amplitudes, which is true for 
T23 as well as for T33, has far reaching consequences on the energy dependence of the 
entire amplitude. Certainly, in order to take this cancelation into account the terms 
corresponding to completely connected diagrams must be included. 

Another general constraint on the form of T23 is provided by the so called Fermi- 
Watson theorem. This requires that the phase of the transition amplitude of the entire 
amplitude is determined by the phases and norms of the rescattering amplitudes in 
the entrance and exit channel and that in the limit of vanishing I SI the overall phase 
of the entire amplitude is equal to that of the elastic scattering amplitude in the exit 
channel. This well established for two-body reactions |13| is easily extended to the 
case of three-body reactions near threshold (see Appendix A) . Explicitly, for the case 
under considerations, this theorem requires that the phase of T 23 must be identical to 
the overall phase of T33 . Both of these constraints are satisfied by the approximation, 
eqn. 9, in a natural way. 

Turning now to estimate the rescattering amplitude, we note first that the kernel 
FF is compact for any complex E and therefore eqns. 11 has a unique solution. 
By solving these equations by iteration , one can calculate the amplitude T33 with 
the required accuracy. To first order T^ l) « [t + Ft}. It is the sum of all 
diagrams shown in fig. 3b. The first three of these are disconnected and represent 
the contributions of the free t. As shown by Amado[jl2| and Cahill[|l4]] taking the 
sum of just these three diagrams would not yield a coherent amplitude and therefore 
would not be a satisfactory solution. The other diagrams of fig. 3b give the first order 
correction due to two subsequent elastic scatterings, they are all completely connected 
and as indicated above play an important role in the coherence mechanism [|T2| , 
The next iteration introduces connected diagrams with three and four rescatterings 
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and so on. Summing the contributions of diagrams up to the corresponding order 
leads, in principio, to T33 with the accuracy required. 

As a first step we restrict the following discussion to FSI corrections due to one 
and two rescattering diagrams only (fig. 3b), i.e. , we calculate to first order 
only. Denoting contributions from double scattering diagrams by Cji it is shown in 
Appendix B that for s-wave production at energies close to threshold, 

C sl = M< fe >(0, 0) Qj t^(0, iL) t j«)( , |L) l iX . (16) 



where Iji are integrals of Kowalski-Noyes |T5|1 half-shell functions over the appropriate 



relative momenta of the interacting pairs. Thus the entire transition amplitude can 
be written as, 

T 2 3~Mg n) (0,0)Z 33 , (17) 
where the FSI correction factor is given by, 

= e *r } (oA + e Qi t ( r ] (o, $-) t ( r n) (0, f-) hi . (is) 

j=l j,l=l l H Z Vl 

A further simplification can be achieved by setting all Iji equal unity, thus neglecting 
off shell effects (see Appendix B). This allows calculating the FSI factor from two- 
body scattering data of the particles involved. To evaluate T33 to second order we 
need to include ( FF ) [t + Ft] terms. The evaluation of such terms is described in 
Appendix C, where it is shown that the effect of the (FF) factor is to scale the first 
order amplitude [t + Ft ] by, 

(FF) « q k t ( r\0,^)<llt < r\0,^) (19) 
w expi( 5k + 5/)sin5feSin5/ . 

Here we have used the standard parameterization, 

'^^mW (20) 

with the 5's being the appropriate s-wave phase shifts. These are related to the s-wave 
scattering length according to, 



1/ tp i? \i ~ • x • x <lk a>k qi a-i / on x 

(FF w sm4sm^ = — — — — . 21 

(1 + q k a k )(l + qi a/) 
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Thus in case of q^ a& q\ a\ <1 , the factor ( FF ) 1 so that the evaluation of T^P 
to first order should be sufficient. 



3 Application 

In the present section we apply the procedure described above to analyze the effects 
of FSI on the energy dependence of tt and rj meson production in NN scattering 
at energies near their respective production thresholds. We demonstrate this first 
for 7r meson production. For the pp — *> ppm reaction, rescatterings may include the 
following sequences : pp — > ppir — > ppir , pp — > pnir + — > ppn° and pp — > dir + — > ppir . 
In considering FSI only the first sequence is included as by definition the block FSI 
includes elastic interactions only. Two-nucleon mechanism contributions dominate 
the pp — > ppr] cross section and any reasonable calculations of the primary amplitude 
must include 7r + production followed by charge exchange. The effects of these other 
channels are assumed to be absorbed in M^ n \ The isosinglet and isotriplet pion- 
nucleon scattering amplitudes are calculated from eqn. 20, using the isosinglet and 
isotriplet pion nucleon scattering length values a\ = 0.245 fm and 03 = — 0.143 fm of 
ref.|16|. The isospin average scattering length of these values is rather small so that 



a priori the overall 7rN interaction effects on the energy dependence are expected to 
be very small. Similarly, the S-wave pp phase shift is obtained from effective range 
expansion using the modified Cini-Fubini-Stanghellini formula JTljj, 

A 



Pi q 

a pp ' 2 ,pp ^ 1 + p 2 k* ' 



C 2 qctnd + 2ri c qh(r, c ) = — + -r pp q 2 - - ^ L * - , (22) 



where a pp and r pp denote the pp scattering length and effective range and, 

e 2 E lab 2 2tt r] c ^ Vc , / 0Q x 

P lab (exp 2Tcr] c - 1) ^ s{s 2 + r, 2 c ) 

The Pi (i=l,2) are functions of a pp and r pp . In what follows we have used the values 
flpp = — 7.82/m and r pp = 2.7 fm of ref. |DJ. To calculate the cross section, the FSI 



factor |Z 33 | 2 , eqn. 18, is multiplied by the invariant phase space and then integrated 
over the appropriate momenta. The primary production amplitude is assumed to be 
constant and is taken outside the integral. The energy dependence of the cross section 
as obtained with the FSI factor calculated to first and second orders are displayed in 
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fig. 4 as a function of Q cm , the energy available in the cm system. The two solutions 
are practically identical. The ( FF ) factor discussed in the previous section, is indeed 
very small so that the sum in eqn. 18 converges and it is safe to compare data with 
the predictions corresponding to first order calculations of IZ33I 2 . Comparison with 
data is made in fig. 5. The predictions with the full FSI correction of eqn. 18 
are shown as the solid line. Those without the irN interaction and with only the 
first disconnected diagrams of fig. 3b are given by the large-dashed and small-dashed 
lines, respectively. The data in the figures are taken from ref.jlj. The curves are 
arbitrarily normalized to the lowest energy data point available. As indicated above 
the s-wave nN interaction is relatively weak so that the FSI correction factor is 
dominated by the strong interaction of the 1 Sq NN system. Our predictions with 
the full FSI correction factor (solid curve) agrees slightly better with data compared 
with those obtained by including the NN force only (dashed line). Nonetheless, even 
in this case, where the overall meson-nucleon interaction is weak, the effects from 
disconnected and completely connected diagrams are comparable. It is our contention 
that the completely connected diagrams of fig. 3b must be included if one wishes to 
interpret the recent precision measurements of 7r production cross section Jl], ||. As 
demonstrated below the importance of these completely connected diagrams becomes 
more apparent in applying our procedure to rj production. 

For 77 meson production there are no direct measurements of elastic r]N scattering 
and the information available are extracted indirectly using ttN — > i]N data. Based 
on ttN and rjN coupled channel analysis around the r)N threshold within an isobar 
model, Bhalerao and Liu|TB[ suggest a value a V N = (0.27 + i0.22) fm for the r/N 
scattering length. With this choice of cl v n the results of our analysis with the FSI 
factor calculated to first order are shown in fig. 6. Here as for pion production, the 
solution with the FSI factor calculated to second order is nearly identical. In contrast 
with 7r production, this process is strongly influenced by the r]N interaction giving rise 
to a sharp enhancement of the cross section very near threshold and thus reproduces, 
rather accurately, the energy dependence of the cross section up to Q cm = 20MeV. 
Note also, that in the analysis presented above, the enhancement factor \Z^\ 2 is 
exactly what is required to keep the primary production amplitude nearly constant 
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over the energy range from threshold to Q cm = 20MeV. 

It is interesting to explore how sensitive the results are for a different choice of 
the a n N scattering length. Since other channels such as the 7T7riV are also opened, 
Wilkin^] has fixed the imaginary part to be ^[a^Tv] = 0.30/m and used the n~p — > rjp 
data directly to determine the real part. This procedure has yielded a scattering 
length a v jsr = (0.55 ± 0.20 + z0.30) fm. Calculations with this value are given in fig. 
7. The energy dependence is mostly sensitive to the real part of a V N and the results 
from our analysis do not support the central value proposed by WilkinQ. It is to 
be noted also, that based on partial wave unitarity relations, an analysis of the 7i~p 
backward differential cross section data yields a value of 39°. 3 ± 4°. 7 for the phase 
of the nN production amplitude [T^] a value in agreement with Bhalerao and Liu[TE| 



but not with Wilkin @]. 

4 Summary And Conclusions 

In summary we have developed an approximate procedure which accounts for final 
state effects in the presence of three strongly interacting particles. In analogy with 
the Watson's treatment |9|] it was shown that, in the the limit of extremely strong 
FSI the amplitude factorizes into a primary production amplitude and an elastic 
scattering amplitude describing a 3 — ► 3 transition. The energy dependence of the 
cross section in the near vicinity of the reaction threshold is determined almost entirely 
by T33 . Based on the coherence formalism of Amado [|l"2]j we argue that the completely 
connected diagrams corresponding to two rescatterings may interfere strongly and 
therefore, play an important role in determining the energy dependence of the cross 
section. In view of the scarcity of knowledge about the meson-nucleon force, the 
approximation is particularly useful near threshold where T33 can be estimated using 
nucleon-nucleon and meson-nucleon scattering data. Because of the strong NN force 
the FSI term dominates. Detailed analysis of tt and rj production via NN collisions 
show that the energy dependence of the cross section can be explained in both cases, 
by taking into account the meson- N interaction also. 

In the present analysis we have not considered ISI corrections as these vary 
slowly with energy and would have little influence on the energy dependence. For 



14 



iViV scattering at bombarding energies of several 100 MeV and above, the nuclear 
and Coulomb distortions are expected to be very small. The Gamov factor which 
accounts for Coulomb distortions is equal unity within 2.5% even for the pion data 
analyzed above. As previously mentioned ISI effects can be included exactly by 
using distorted waves for the incoming two-proton state. 

Finally, the primary production amplitude is not treated explicitly so that the 
cross section scale was arbitrary. This deficiency can be removed by incorporating 
in eqn. 17 a primary production amplitude from one of the existing perturbative 
models^], || [Tj. Such model calculations were performed for 77 production, taking the 
??N interaction in the primary amplitude to be consistent with that occurring in the 
FSI . They are found to reproduce the energy dependence as well as the scale of the 
cross section and will be published elsewhere [|20|] . 

In view of the large effects due to the rjN force the analysis presented above could 
well be used to study the rjN force itself for which no direct experiments are possible. 
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Appendix A 



The Fermi- Watson theorem [|T^] can be extended to three-body reactions by ap- 



plying unitarity. Let T 23 denotes the entire amplitude and let t^) and denote 
elastic 2 — > 2 and 3 — ► 3 amplitudes, respectively Then by unitarity, 

3{T 23 } = -7T J2 if 5 ( E - E *) T k - * E T ™ 5 i E - E 3 ,)T^ . (24) 

2' 3' 

The 5 functions force the two body and three body elastic scattering amplitudes, t^) 
and T3/3 to be on mass-shell. 

In angular momentum representation the amplitudes of eqn. 24 can be written in 
the form, 



Tij = \Tij\ expi6ij . (25) 

Note that 8^2 stands for two-body phase shift but #3/3 is only related to two-body 
phase shifts of the three interacting pairs and in itself is not a measurable quantity. 
Then substituting these into the unitarity condition, eqn. 24, leads to the following 
constraint on the <5's, 



S{|T 23 | expi5 23 } = -71~E 1^22' 1 1^2'3 I eX P *(<^22' - S 2 ' 3 ) 

2> 

l T 23' n) I l T 3'3 I eX P «(<*23' - 5 3 ' 3 ) • (26) 

3' 

In the case of a single channel with the angular momenta of the final pairs being 
zero, the phase £ 23 from eqn. 26 is, 



Sin S 23 = |*22 I sin ^22 - 1^33 I sin <W 

[|4r } | 2 + \Tt ] \ 2 - 2|42 n) ||T 3 ( 3 n) | cos(5 33 - 5 22 )] 1/2 • (27) 

Thus the phases of the rescattering amplitudes and their norms, determine the 
overall phase of the 2^3 transition amplitude, just as the case is for a two body 
reactions. In the limit of very weak I SI (vanishing tj^P), 
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sin 5 23 « sin 5 33 , (28) 

so that the overall phase of T 2 3 becomes identical to that of the elastic scattering 
amplitude Furthermore, the amplitudes and phases of the strong rescatterings 

in the entrance and exit channels are related through eqn. 26 by, 

(sinfe + tt |4 e 2 °|) = |T 3 ( 3 e °| (sinfe + vr |T 3 ( 3 e °|) . (29) 

For a real production process the amplitudes t%P and T33 must be independent 
so that the lhs and rhs of eqn. 27 must be equal to a constant. In the limit of very 
weak I SI this constant must vanish so that, 

\T 3 f\ =sin<W7r • (30) 

Formally, this is identical in form to the expression assumed by Watson for two- 
body processes. 

Appendix B 

We evaluate the contribution of diagrams with two rescatterings. As an example 
consider the diagram shown in fig. 8 (crossed lines denote on mass-shell states). In 
the CM system, in each of the three-body states of the diagram, there are only two 
independent momenta which are taken to be pj the momentum of particle j and qj 
the relative momentum of the remaining / — k pair. Adopting the convention that 
(j, /, k) are always cyclic, q., is defined as 

= (m k pi - mip k ) 
qj ( mi + m k ) 

In terms of these momenta the total kinetic energy is , 

Hq= pL + A- . (32) 
2//j 2// Zfc 

Here /i^ is the reduced mass of the / and k particles and fj,j is the reduced mass of 
the j particle, i. e. , 

mim k rrij^mi + mk) 

(mi + m k ) (rrij + m t + m k ) 
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Now the contribution of the diagram, fig. 8, can be written in the form, 

C31 = /W(E-^-^)MW( P ;,q;) (34) 



n" 2 

dp'^q'^p's " p£)fs(<& OS; ^ + " fM 

^3 

//2 //2 2 

[e + - f - - 77^ ]~ 1( Kpi - Pi))ti(q , i', qi; £ + ' 1 



2/'i -/'i'.'. -/': 
The integration over p" and p' 3 are immediate and we obtain, 

r r>" 2 a' 2 

C 31 = d^8(E-^-^)M^(pl^) (35) 



[ 2/i 23 J U41 ' 41 ' 2/V 

Here £3 is a half off mass-shell two-body matrix element with only the q 3 ' being off- 
shell momentum. (From a view point of invariant perturbation theory only one leg 
of tj is off mass-shell so that it behaves like a particle form-factor). Expanding tj in 
partial waves and introducing the Kowalski-Noyes half-shell function[|15[ fj(X, q,, q") 
allows writing these for a partial wave A as, 

tj{\ q„ <# E+) = tf n) (X, fj(X, q,, t$) . (36) 

Then for slowly varying and at energies near threshold one obtains, 

Csi = M^(0, 0) g 3 4° n) (0, ^-) *i° n) (0, hi , (37) 

^12 ^23 

where J31 is an integral over the Kowalski-Noyes half-shell functions, 

1 2 //2 

^31 = 7^2 fdtfMO^&lE-^-^-'fifrqx,*® , (38) 

(2tt) 2 /x 2 3 93 J 2/ii 2/i 23 



where 



„ mx „ m 2 (m 1 + m 2 + m 3 ) 

q 3 = — qi - 7 — - — T7 — - — v Pl • ( 39 ) 

mi + m.2 \m\ + m2){mi + m 3 ) 
If we neglect off mass-shell effects the integrals hi reduce to unity. 
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Appendix C 



We consider here second order contributions to the transition amplitude. As 
indicated in the text these are obtained by multiplying first order contributions with 
a factor (FF). We attempt demonstrating this for the diagram shown in Fig. 9. This 
can be written as, 

D 3 i3 = / dp' 3 dc^dp'l dq? dpg dqg S(E - |^ - (40) 

M^(p 3 , q' 3 ) t 3 (&, <& E+ - £-) 5(p> 3 - p 3 ') 

^3 

//2 „//2 ///2 

[E+ - EL. - fc^l-i <5( p " _ p '") tl ( q ", q "' ; £+ - 

///2 ///2 2 

[£ + " " l^]- 1 W - Pa) *3(qS', q 3 ; £ + - ■ 
Integrating over p 3 , p" and p 3 ' and arranging terms leads to, 

„///2 „2 1 „ 1 

£313 = { ?i *i° n) (0 , |M qs tt\0 *)J- dq« ' ( 41) 

2/i 2 3 2//12 gi g 3 •/ LB+ - ^ ^j— 

f, ,,, /i(0,q?, 0/1(0,^,03) X 

111 12 12 

[ rfq 3 MW(p; q 3 ) - |?_ _ J3_ )t3(0? ^_) /3 ( , q 3 , q») . 

We now notice that the last integral is the first order contribution from the discon- 
nected diagram (9. a) with a single scattering block between particles 1 and 2 while the 
expression in the curly bracket is just the (FF) factor. An order of magnitude of this 
factor can be obtained using on mass-shell values for the Kowalski-Noyes functions. 
This leads to, 

W2 2 

(FF)*> qi t[° n \0 , |H q 3 i on \0 , -*-) . (42) 

^23 ^12 

Similar expressions can be calculated for the completely connected diagrams (3.b) as 
well. 
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Figure 1: Graphical representation of the transition amplitude. The block 
represents the primary production mechanism of the reaction and contains all possible 
inelastic transitions from a two-particle state to a three-particle state (see text). The 
ISI and FSI blocks represent elastic rescattering in the entrance and exit channels, 
respectively. 
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(a) (b) (c) 



Figure 2: Elementary inelastic(2a) and elastic(2b) interactions. A wavy line repre- 
sents the primary interaction, such that if it were zero the reaction ab — > cd would not 
occur. A dotted line represents an elastic interaction. The disconnected diagram(2c) 
describes non-interacting particles. 
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(b) 



Figure 3: Graphical representation of elastic amplitudes: a) the amplitude t^2 f° r 
ab — > ab, b) the amplitude for cde — > cde; the blocks represent two-body scat- 
tering amplitudes in three-body space similar to the one defined in 3a. 
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Figure 4: The pp — > ppn° cross section vs. the energy available in the cm system, 
Qcm- The solid line represents predictions with the FSI correction factor calculated 
to first (T33 ~ [t + Ft]). The dashed line is that obtained with the FSI taken to 
second order,ie., T 3 ( |° w (1 + FF)[t + Ft] (see text). 
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Figure 5: The ir meson production cross section vs. the energy available in the cm 
system, Q crn . The predictions with the full FSI correction factor of eqns. 20 are 
shown as the solid line. Those without the meson-nucleon interaction and with the 
disconnected diagrams of fig. 3b are given by the large-dashed and small-dashed lines, 
respectively. All curves are normalized to the lowest data point. 
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Figure 6: The rj meson production cross section vs. Q cm - The r/N scattering len 
is taken from ref. [18]. (see captions of fig. 5) 
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Figure 7: The rj meson production cross section vs. Q cm - The rjN scattering length 
is taken from ref. [3] (see captions of fig. 5) 
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Figure 8: Typical completely connected diagram with two subsequent rescatterings. 
Crossed lines denote on mass-shell states. 
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Figure 9: An example of diagrams contributing to the T^P amplitude to second 
order : (a) a disconnected diagram which contributes to first order, (b) a completely 
connected diagram obtained from diagram (a) with the FF factor (see text). 
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